Evolution of magnetic oxygen states in Sr-doped LaCoC>3 
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Magnetism in Lai-^Sr^CoOs as a function of doping is investigated with X-ray absorption spec- 
troscopy (XAS) and X-ray magnetic circular dicrhoism (XMCD) at the O K-edge, and corresponding 
first principles electronic structure calculations. For small x, the spectra are consistent with the 
formation of ferromagnetic clusters occurring within a non-magnetic insulating matrix. Sr-induced, 
magnetic O-hole states form just above and grow with increasing Sr doping. Density functional 
calculations for x=0 yield a non-magnetic ground state with the observed rhombohedral distortion 
and indicates that doping introduces holes at the Fermi level in magnetic states with significant 
O 2p and Co t2 9 character for the undistorted pseudo-cubic structure. Supercell calculations show 
stronger magnetism on oxygen atoms having more Sr neighbors. 



Doping in complex oxides is a fundamental route for 
control of the many degrees of freedom. However, de- 
pending on the nature of the underlying ground-state, the 
impact on the local electronic structure depends strongly 
on the balance between on-site energies and the bonding 
to oxygen P, For example, doping in the Mott regime 
directly changes the valence state of the transition metal 
ion (e.g. manganites) while doping in the extreme of the 
charge transfer regime (e.g. cuprates) results in forma- 
tion of holes on the oxygen lattice. The case of cobaltites 
lies in the cross-over region between these two extremes 
and a better understanding of this regime may finally 
answer the decades-old question of why magnetism in 
Lai^zSr^CoOs (LSCO) is so unusual. 

LaCoC>3, is a non-magnetic, small gap semiconductor 
at low T, but at increased temperature and particularly 
with Sr doping on the La sites, the material exhibits a 
wealth of unusual magnetic properties associated with Co 
sites acquiring a substantial magnetic moment (^1.3/is). 
For x > 0.18 the crystal becomes ferromagnetic (FM) and 
metallic, 0-Q — see Fig. 5 of Ref. for a phase diagram. 
With increasing x, the rhombohedral structural distor- 
tion gradually decreases and the crystal approaches the 
cubic perovskite structure Q. 

Most studies of LSCO have used localized models 
which assume that the dominant interactions are the 
crystal field splitting and the exchange interactions E ea; ; 
these interactions split the Co 3d energy levels into nar- 
row t2 g and e g multiplcts and also split the spin up and 
spin down states. Within such models there is exten- 
sive debate about the spin states, with low spin (S=0, 
t| g e g)@ at low T and a possible mixture of low spin 
and either an intermediate spin (S=l, t^eglfiot or a 
high spin state (S=2, t| g e 2 ) at higher T or with Sr dop- 
ing. However, we argue that this localized picture is not 
appropriate — the e g states interact with O 2p states and 
are spread over a large energy range (10 eV), from low 
energy bonding states to high energy anti-bonding states. 



In Co K-edge X-ray absorption near edge spectroscopy 
(XANES) studies^ E3 the Co K-edge shifts < 0.15 
eV as x increases to 0.3, in stark contrast to other Co 
systems and many manganites, for which the edge shift 
is roughly 3 eV per valence unit(l3- 15|. Further, both 
diffraction studies fl6L [Trj] and extended X-ray absorption 



fine structure (EXAFS) results[12j show that the Co-0 
bond length is nearly independent of x for LSCO; con- 
sequently, the bond-valence sum model also indicates no 
change in Co valence — i.e. the 3d electron configuration 
remains close to d 6 . This raises the question: where do 
the Sr-induced holes go? We argued that a large fraction 
of the holes are on the O atoms [12J , which can be directly 
probed at the oxygen K edge. Previous XAS studies at 
the O K-edge show an evolution with hole doping [lj| Gil, 
but existing X-ray magnetic circular dichroism (XMCD) 
studies arc limited. I2C 
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Thus, to date there has been 
no systematic study of the low hole doping regime to fol- 
low the connection between local structure and changes 
in electronic and magnetic states. 

Here, we present a detailed doping-dependent study 
of the oxygen hole states in LSCO single crystals, using 
polarization-dependent XAS to track the first O pre-edge 
peak (Co-0 ligand hole states). It is clear that a mag- 
netic state with a large O component, forms with doping 
and the increase in O absorption is tied to a large frac- 
tion of the doped holes residing on the oxygen site. The 
strong XMCD shows directly that there is a non-zero or- 
bital moment on oxygen (and correspondingly, a net O 
magnetic moment). The analysis shows the moment on 
oxygen is parallel to the Co moment and increases in 
magnitude with doping. To understand these results and 
determine the magnetic structure of O, complementary 
theoretical studies have been carried out; they show that 
the 3d electron density of states (DOS) are spread over 
a large energy range, with significant e g character for 
bonding states 6 eV below Ep. They also show that the 
(empty) O DOS just above the Fermi level increases with 
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FIG. 1: XAS at the O K-edge. The XAS data were normalized 
at 550 eV and self-absorption corrections were applied. The 
data agree well, except at the beginning of the edge (from 527- 
529 eV) where, with increasing x, a large increase in the XAS 
is observed, corresponding to the addition of empty states. 
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FIG. 2: Difference in the XAS (top) and XMCD (bottom) 
at the O K-edge. A significant XMCD signal is only ob- 
served over this same small energy range (shaded region) 
where change in the XAS occurs and XAS(x)-XAS(0) is pos- 
itive. Above 529 eV, the differences become slightly negative. 

x, and that the O atoms develop a moment, with both 
spin and orbital components. Most importantly, the O 
atoms are no longer equivalent; their properties depend 
on the number of nearest Sr neighbors. 

XMCD and XAS data for the O K-edges were mea- 
sured at the Advanced Photon Source beamline 4-ID-C. 
A superconducting magnet provided magnetic fields up 
to ±5T. The energy resolution was ~ 0.15 eV. Small sin- 
gle crystals were mounted on an electrode using Ag paste 
and scraped with a diamond file to remove surface con- 
tamination. The X-ray beam was oriented at 45° to the 
sample surface and the fluorescence detector was at 90° 
to the beam. 

The XAS fluorescence data at the O K-edge are shown 
in Fig. [TJ A small constant background was subtracted 
from the XAS and then the data were normalized to 1 at 
energy 550 eV. Finally, the data were corrected for self 
absorption following the method in the computer pro- 
gram FLUO[i2|. The main results for the XAS data as 
x increases are a shift of the leading part of the O K- 
edge to lower energy and an increase in absorption. As 
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FIG. 3: Integrals of the difference in O XAS (Fig. [2]) from 526- 
529 eV (magnetic holes) and from 529-532.5 eV (non-magnetic 
holes) as a function of x; note that the change is not linear. 
For x < 0.15, the new magnetic holes are associated with FM 
clusters; over the range 0.15-0.22 where the material becomes 
fully FM, higher energy holes with no orbital momentum are 
converted to lower energy magnetic states. 

has been noted previously [23| and in our theory results 
below, the weak core-hole interaction at the O K-edge 
means that the XAS closely follows the shape of the O- 
projected density of states (DOS). This allows us to make 
a direct connection with the evolution of electronic struc- 
ture. 

The magnetic state on O is probed by XMCD 
where fi + and fi~ arc the absorption coefficients for ± cir- 
cularly polarized X-rays with the sample in a B-ficld high 
enough that the magnetization is close to saturation) as 
shown in the lower part of Fig. [5J A significant signal 
is only observed over a short energy range (~ 2 eV) at 
the beginning of the O K-ed ge, indicating a non-zero or- 
bital moment on the O site 24j . The amplitude of the 
XMCD signal increases with x and is surprisingly large 
above x = 0.19 compared to measurements of O XMCD 
in other systems; the formation of the FM metallic state 
appears to arise from the decreasing rhombohedral dis- 
tortion with x. However, the fact that only a small part 
of the O XAS shows a magnetic signal indicates that not 
all the O hole-states are magnetic. 

To visualize how the electronic structure evolves with 
x, the Sr-induced change in XAS defined as XAS(x) - 
XAS(.t = 0) is shown in the top of Fig. [5] Sr induces 
a strong increase in the absorption at lower energy due 
to the formation of additional unoccupied oxygen ligand 
hole states. For low x, this new XAS peak shifts to lower 
energies and increases in amplitude as x increases, with 
no other changes in the XAS up to <~ 533 eV. For x > 
0.15, there is no further shift, but the amplitude contin- 
ues to increase with x and the XAS difference becomes 
negative above 529 eV {i.e. the O-DOS decreases from 
529-532 eV above the metal-insulator transition (MIT)). 
The shaded regions in Figs. [1] and [5] illustrate that the en- 
ergy range for which the XMCD is observed corresponds 
to the same energy range over which the additional peak 
in the XAS develops. 

This increase in empty O-DOS is consistent with Sr 
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FIG. 4: 2x2x2 supercell for a Sr doping of 3/8. Sr (large, 
light blue); La (large, dark green); Co (medium, blue). Four 
types of (small) O atoms with varying Sr neighbors: (red, 
left slant); 1 (magenta, vertical); 2 (pink, horizontal); 3 (yel- 
low, right slant). 
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doping adding holes to the hybridized O bands. Note, 
however, that the density of the oxygen states with 
quenched angular momentum (i.e. above 529 eV) de- 
creases with x for x above ~ 0.18. The induced O mag- 
netic holes (DOS integral from 526-529 eV) and non- 
magnetic holes (from 529-532 eV) are shown in in Fig. 
[3l There is a clear change in these integrals near x = 
0.18 where the MIT occurs and the system changes from 
a spin glass to a FM state. 

To understand this evolution, we utilize doping- 
dependent electronic structure calculations. Recent stud- 
ies suggest that a first principles density functional ap- 
proach from an itinerant electron viewpoint can pro- 
vide insights about the behavior of LSCO.jH-Ell This 
is a homogeneous approach, which is typically applicable 
for itinerant or free-electron systems. We use the elec- 
tronic states obtained from first principles band structure 
calculations to calculate theoretical XAS and XMCD 
spectra. For the electronic structure, we use a full- 
potential linearized augmented plane-wave (FPLAPW) 
method [11] with a local density functional. [29l l3fjj We 
used R M T^max = 8.0 and R M t = 2.2, 2.2, 1.9 and 1.6 
a.u. for La, Sr, Co and O respectively. 

Because of the strong short range electronic interac- 
tions, the virtual crystal approximation (VCA) is not 
appropriate. Consequently, we use a 2x2x2 supercell 
formed of eight pseudo-cubic unit cells, with Sr concen- 
trations 0, 0.125, 0.25 and 0.375 as shown in Fig.g] We 
initially kept the cubic structure even for the lower con- 
centrations, since preliminary calculations indicated that 
the spectral features were more sensitive to the local 
atomic arrangement compared to the overall cell geom- 
etry. For the supercells, 40 k-points were used for the 
iterations to self-consistency. 

The strong interactions between O 2p and Co 3d (par- 
ticularly c g ) orbitals form well-separated bonding (-6 
eV) and antibonding (-1 to 4 eV) states — see the or- 
bital decomposed DOS in Fig. [3] A local model with e g 



FIG. 5: Orbitally decomposed DOS for Co 3d (t2 9 and e 3 ) 
and O 2p states from LDA calculations. For each panel top 
is spin up, bottom, spin down. The shaded region for E>Ef 
is expanded on right. 



mainly above the t2 g states is inconsistent with these 
DOS. The generalized gradient approximation (CCA) 
and local density approximation (LDA) density function- 
als yield very similar DOS and only the LDA results are 
plotted. Our calculations for undoped LaCoOa are essen- 
tially in agreement with those of Ravindran et al. 26|, 27 1 . 
However, although the GGA[3(| is more sensitive to the 
electronic correlations and yields a nonmagnetic ground 
state with experimental rhombohedral distortion (Co-O- 
Co bonding angle is 163.4°) LDA needs stronger distor- 
tion (Co-O-Co bonding angle is 158.4°) to yield the non- 
magnetic state. The transition from the nonmagnetic to 
magnetic ground state (with ~4.3/is moments on Co) is 
not unlike the traditional local transition state picture, in 
the sense that about one Co d electron is promoted from 
spin down to spin up with very little cost in energy, and 
is very sensitive to the rhombohedral distortion. Using 
the experimental value for the rhombohedral distortion, 
the CCA calculation shows a ~ 3.2 meV/Co total energy 
difference between non-magnetic and magnetic states. 

We also extracted the theoretical O K-edge (mainly a 
measure of the oxygen partial DOS (PDOS)) and the O 
XMCD signal, both as a function of x, as seen in Fig. 
El These theoretical results agree quite well with the 
experimental results shown in Figs. Q] and [2J though it is 
important to keep in mind that the energies are relative 
to Ep, so the XAS does not shift to lower energy as 
is observed experimentally, but the increases are quite 
similar and, again, we see there is only a large XMCD 
signal over the same energies for which the XAS has a 
large increase. The lack of an amplitude change with 
x is likely because the rhombohedral distortion was not 
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FIG. 6: Theoretical O K-edge (top), showing a significant 
increase at the beginning of the edge, and calculated O PC- 
edge XMCD signal (bottom). 
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FIG. 7: Calculated hole fractions: magnetic O-hole fraction 
(red square), non-magnetic O-hole fraction (green triangle), 
total O-holes (blue circle), total Co-holes (gray star). 

included for small x. 

Further important details are obtained by calculating 
the Sr-induccd hole fraction on Co and O atoms, using 
similar energy ranges as for Fig. |3l In Fig. we plot the 
induced hole density as a function of x; the calculated O 
hole density is comparable to but slightly larger than that 
on the Co atoms. Thus, surprisingly the induced holes 
are nearly equally distributed over the Co and O atoms, 
in good agreement with the experimental findings. 

The calculations of the individual O magnetic moments 
are even more unusual and show that the O atoms in the 
supercell are not equivalent. Their properties depend 
on the number of neighboring Sr atoms, which can vary 
within the supercell (from to up to 3 for £=0.375): the 
larger the number of Sr neighbors, the larger the calcu- 
lated magnetic moment (from 0.067/xb for O with Sr 
neighbors to 0.121/xs for those with 3). Thus, these cal- 
culations show the limitations of the VCA; the system is 
not homogeneous at the local level. A similar situation 
occurs for Co, showing different Co moments on differ- 
ent sites, when calculated using a trigonal distortion of 
the basic unit cell, since the cubic 2x2x2 supercell is 
too small to discriminate between different Co. Another 
interesting feature from the calculated O-PDOS (Fig. [5]) 
is that it provides a simple explanation as to why there 
are magnetic and non-magnetic holes; just above Ej? the 



spin up and spin down O-PDOS are not equal so those 
energy states will have a net moment. At higher energies 
the spin up and down O-PDOS are nearly equal and will 
have little moment . 

In summary, we find that it is better to view these 
cobaltites in terms of bonding and anti-bonding states 
rather than split t2 9 and e ff states. Some of the e g sym- 
metry states form bonding states and are filled; conse- 
quently, there is never the situation where some e ff are 
empty A mixture of filled e g and t2 9 states may explain 
the wide range of conflicting experimental results about 
the Co spin state. Further, the broad range of the e g DOS 
precludes a Jahn- Teller (J-T) interaction from being im- 
portant. This likely explains the lack of a significant J-T 
distortion for the Co-0 bonds. 
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The experimental changes in the O XAS and XMCD 
are not linear with x and appear to be correlated with 
the overall magnetism in the system, which changes from 
glass-like with FM clusters below 18% to metallic FM 
above 18%. Calculations qualitatively agree quite well 
with both the XAS and XMCD at the O K-edge, and 
indeed find a small but significant spin magnetic moment 
on O atoms. Further, the theoretical results show that 
the properties of the O atoms are not all equivalent, and 
have different moments; confirming that the VCA is not 
adequate for describing systems such as LSCO, for which 
the electron density and oxygen moments are sensitive 
to nearby Sr positions. With part of the doped holes 
being on the O atoms, the Co-0 Coulomb attraction that 
usually contributes to a shorter metal-0 bond when holes 
are added, is reduced and likely is a major reason as to 
why there is so little change in the Co-0 bond length 
with doping. 

Note: a recent explosive high field measurement up 
to 500 T on single crystal LaCo03 indicated a magnetic 
state appearing above 140 T with a magnetic moment of 
~ 1.4 fiB, which could correspond to the excited state 
with the same moment obtained by our band structure 
calculations. j3lj 

The XAS and XMCD measurements were carried out 
at the Advanced Photon Source at Argonne National 
Laboratory and is supported by the U.S. Department of 
Energy, Office of Science under grant No. DE-AC02- 
06CH11357. Work at the Ames Laboratory was sup- 
ported by the US Department of Energy, Basic Energy 
Sciences under grant No. DE-AC02-07CH11358. 
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